home *** CD-ROM | disk | FTP | other *** search
/ Info-Mac 4 / Info_Mac IV CD-ROM (Pacific HiTech Inc.)(August 1994).iso / Science / RLaB / examples / primes.r < prev    next >
Text File  |  1994-04-25  |  340b  |  23 lines
  1. //
  2. // An example that finds all primes less than limit
  3. //
  4.  
  5. primes = function (limit) {
  6.   local(prime, cnt, i, j, k);
  7.  
  8.   i = 1; j = 0; cnt = 0;
  9.   for(k in 2:limit) 
  10.   {
  11.     j = 2;
  12.     while(mod(k,j) != 0) {
  13.       j++;
  14.     }
  15.     if(j == k) {            // Found prime
  16.       cnt++;
  17.       prime[1;i] = k;
  18.       i++;
  19.     }
  20.   }
  21.   return prime;
  22. }
  23.